Transferring entanglement to the steady-state of flying qubits 
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The transfer of entanglement from optical fields to qubits provides a viable approach to entangling 
remote qubits in a quantum network. In cavity quantum electrodynamics, the scheme relies on the 
interaction between a photonic resource and two stationary intra-cavity atomic qubits. However, it 
might be hard in practice to trap two atoms simultaneously and synchronize their coupling to the 
cavities. To address this point, we propose and study entanglement transfer from cavities driven 
by an entangled external field to controlled flying qubits. We consider two exemplary non-Gaussian 
driving fields: NOON and entangled coherent states. We show that in the limit of long coherence 
time of the cavity fields, when the dynamics is approximately unitary, entanglement is transferred 
from the driving field to two atomic qubits that cross the cavities. On the other hand, a dissipation- 
dominated dynamics leads to very weakly quantum-correlated atomic systems, as witnessed by 
vanishing quantum discord. 

PACS numbers: 03.67.Mn, 03.65.Ud, 42.50.Pq 



I. INTRODUCTION 

Entanglement is an invaluable resource for many im- 
portant tasks in quantum information processing l5|j , 
including metrological purposes d, 0] ■ Schemes for the 
generation and distribution of entanglement have been 
designed and implemented, in the past 20 years, in a 
number of physical systems. The quality of these strate- 
gies has considerably increased over recent years, making 
the creation of multipartite entangled states of a few el- 
ements a reality pfrjlQj and paving the way to the near- 
future realization of networks of distributed quantum 
nodes for quantum communication and computing 
The main challenge, in this context, is the achievement 
of reliable interfaces between information carriers having 
different natures. 

This problem has long been investigated, both theoret- 
ically and experimentally, and various solutions for the 
achievement of controllable interactions between static 
local processors and flying information carriers have been 
designed [5J. Among them, a promising one for its tech- 
nologically realistic nature and its flexibility is embodied 
by the transfer of quantum correlations from light to mat- 
terlike systems [llT[l7| . Broadly speaking, this paradigm 
for entanglement distribution via light-matter interfaces 
requires the availability of entangled-light resources and 
the ability to perform local light-matter interactions that 
pass the quantum correlations (or part of them) to ini- 
tially separable local matterlike systems. Originally de- 
vised for cavity and circuit QED settings |hO - [l3| . this 
approach has recently been extended to mechanical sys- 
tems interfaced to light [l8j and quantum many-body 
systems [l9| . 

In most of the above-mentioned schemes based on 
cavity-QED technology, it has been assumed that multi- 



level atoms are located at a fixed point within a cavity, 
where they are coupled to an antinode of the cavity-field 
standing wave. In the actual situation, the atoms are 
either flying through a cavity [2(| or trapped within an 
intra-cavity optical dipole trap, yet still moving within 
it [2l[. Current records for atomic optical traps reach 
trapping times as long as 30 s [22} . However, the man- 
agement of a network of single atoms trapped in distant 
cavities embodies a considerable challenge. Moreover, 
although efficient ways to switch the cavity-atom inter- 
action exist, such a configuration makes it hard to arrest 
the evolution of the inter-atom entanglement so as to 
achieve a desired value. 

In this paper, we propose a simple strategy to bypass 
both such difficulties. We study entanglement transfer 
to two flying atomic qubits crossing two remote cavities 
that are driven by an entangled resource. The scheme 
relics on the high level of synchronization that can be 
arranged for the passage of two different atoms across 
remote cavities. Moreover, as the light-atom interaction 
is turned off after the atoms exit the cavities, it is possible 
to arrange for steady-state atomic entanglement. 

As entangled resources, we consider both NOON states 
and entangled coherent states [23] (ECSs), which are 
important non-Gaussian resources of experimental rele- 
vance (five-photon NOON states have been produced re- 
cently [HI, while ECSs with an amplitude of about one 
photon can be produced, using a beam splitter, from the 
coherent superpositions of coherent states described in 
Ref. 25]). Both classes of states have broad prospec- 
tive applications in either quantum metrology [261 ] or 
fundamental research [27j ■ We study the efficiency of 
the entanglement-transfer scheme in various dynamical 
regimes, highlighting the effectiveness of steady-state en- 
tanglement distribution in the quasi-coherent case corre- 
sponding to the good-cavity limit, while stating a no-go 
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result for the regime dominated by cavity dissipation. 

The remainder of the paper is structured as follows. In 
Sec. |n] we describe the system we address and describe 
the preparation of the entangled cavity fields. Section ITO1 
is devoted to the analysis of the relative performance of 
the scheme under the two classes of entangled resources. 
In Sec. HVI we address the dissipation-dominated dynam- 
ics, showing the impossibility of entanglement transfer. 
Finally, we draw our conclusions in Sec. [Vj 



II. THE MODEL 

Let us introduce here the scheme that we address in 
this work, which is shown in Fig. [T] Two freely propa- 
gating field modes, labeled in the following a and 6, drive 
two remote single-mode cavities, which we call A and B, 
respectively. The state of the driving fields will be speci- 
fied, when needed, later on. The cavities are prepared in 
their vacuum state. The coupling between each external 
driving field and the corresponding cavity is modeled as 
a beam splitter (BS) with transmittivity (reflectivity) T 
(R = 1 — T) [28| that is quantitatively determined by the 
quality factor of each cavity [13j] . The state of the cavity 
modes that result from the driving process is given by 



(1) 



where B Aa and B B b are the BS operators. As a re- 
sult, in general, the cavity fields become entangled. Two 
two-level atoms (qubits) with energy eigenstates |0), and 
|l)j (i=l, 2) and Bohr frequencies equal to the frequency 
of the cavity fields pass through their respective cavi- 
ties. Here we consider two cases: (a) the simultane- 
ous free- falling of the atoms across the cavity (i.e., the 
qubits are accelerated by gravity) [2(| ; (b) the constant- 
speed passage of the atoms, which can be realized exper- 
imentally by embedding an optical lattice loaded with 
a single atom into the cavity [22|, H^] or using opti- 
cal traps for atomic confinement and transport. For 
short time intervals, the coupling strength between each 
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FIG. 1: (Color online) Scheme of the protocol: Two remote 
and identical cavities are driven by a two-mode quantum cor- 
related field, which is coupled to each cavity via a leaky mir- 
ror. Two qubits are then sent to pass simultaneously through 
the cavities along the vertical (a;) direction from the tops of 
the geometrical centers. 



atom and the corresponding cavity field can be consid- 
ered as a constant, and the local qubit-field interaction 
can be treated through the standard (resonant) Jaynes- 
Cummings (JC) model that reads, in the interaction pic- 
ture with respect to the free energy of qubits and cavi- 
ties, as ^fc, (t)=«ni(x(t) , y (*)) (fci | 1>j <0| +*J 1 0>j <1] ) (with 
j=l, 2 and ki=A, k2=B). Here, f2(x(i), y(t)) is the qubit- 
field coupling frequency, which depends on the position 
of an atom within the respective cavity and the form of 
the field mode. As a simplifying assumption that does 
not affect the generality of our analysis, we consider the 
atoms falling perfectly vertically. This means assuming 
y=0 in our model (we take the origin of each reference 
frame at the geometrical center of the respective cavity). 
For a transverse Hermite- Gauss mode, Q(x(t),0) takes 
the form Si 



Q(x(t),o) = n 



|W M ,„(a;(t),0)| 
|*0,o(0,0)| ' 



(2) 



where £! is a constant and ^ u>v (x(t), 0) is the mode func- 
tion of the cavity, 

* 2 w f \/2x(t)\ 
¥„, (a:(t),Q) = C u , v e H u I 2— ^ J . (3) 



-l 

'u,v 



Here, we have introduced the coefficient C u 
(2 u 2"u!w!) 1 /2( w 2 7r / 2 ) 1 /2 ; w hile H p (x) is the Hermite 
polynomial of order p and argument x. The waist ujq of 
each field mode is determined by the radius of curvature 
of the mirrors and the cavity length poj . 

Let us first assess the unitary dynamics resulting from 
neglecting any source of noise and losses. The associated 
unitary operator is Ui2AB{t) = U\A{i)®U2B{t), where 
Ujkj(t) = exp[— i J Q Hjkjit^dt' /h}. The dynamics of the 
two flying qubits is, however, nonunitary and described 
by the reduced density matrix 

Pia(<) = ^ AB [U 12AB (t) Pl2 (0) ® PAB (0)ul 2AB (t)}. (4) 

In order to fix the ideas, we assume that the qubits are 
initialized in 1 00) 12 (00 1, which is usually the prepara- 
tion for which the entanglement-transfer process is op- 
timized [3l[ . 

As for the entangled resource, in this paper we will 
concentrate on two families of non-Gaussian states of 
experimental significance, namely, the NOON and ECS 
families. The first family is described as 

|NOON) a6 = -i=(|N0) ah + |0N) afc ) (NGN), (5) 

while the ECS representative that will be used in our 
study is 



|ECS) Qb =A4 e (|aO) ab + |0< fc ), 



(6) 



where Af a = [2(I+e~' Q )]~ 1//2 is the state normaliza- 
tion. In what follows, we study quantitatively the 
entanglement-transfer scheme performed using each of 
such states. 
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III. PERFORMANCE OF THE 
ENTANGLEMENT-TRANSFER PROCESS 

Let us begin with the NOON family, so that the initial 
state of the pump-cavity system is |NOON) a f,iX)|00)AB- 
The field that is thus prepared within the cavities is 32 1 



1 N 



P AB =2 ^ [ m jR m T m (\m0) AB (m0\ + \0m) AB (0m\) 

m=0 



T 

— (\N0) AB (0N\ + \0N)ab(N0\). 



(7) 



An interesting point that will be reprised and explained 
later on in this paper is that for N > 1 the off-diagonal 
elements of the two-qubit density matrix are all exactly 
null, so that no entanglement is transferred. On the other 
hand, for N = 1, the evolved atomic state takes the form 
(in the ordered basis {|11), |10), |01), |00)} 12 ) 



P12W 



/o 

TS/2 TS/2 

TS/2 TS/2 

\0 l-TSj 



(8) 



where S = sm 2 [p{t)} and fj,(t) = £ il(x(t'),0)dt' . It 
should be noted that the impossibility of entanglement 
transfer for N > 1 is not related to the fact that a 
single-photon entangled state violates a Clauser-Horne 
Bell inequality more than any other NOON state with 
N > 1 [331 ] but to the intrinsic structure of correlations 
in the resources at hand and the entanglement-transfer 
process ruled by the JC model. 

We now consider another state of the driving field, 
namely, an entangled coherent state [23j , which is known 
to offer advantages over the NOON class. For instance, 
in some cases ECSs show a remarkably improved sensi- 
tivity for phase estimation as compared to that of the 
NOON states [34J. Such superiority, as we will see, does 
not extend to the entanglement-transfer paradigm per- 
formed using the JC model for the local atom-field in- 
teraction, as ECSs exhibit lower entangling power than 
NOON states. 

An observation that allows for the understanding of 
the depleted entangling capacity of ECSs comes from re- 
alizing that the latter can be seen as a superposition of 
NOON states as 0, H| 

00 n 

|ECS) Qb = A/- Q e~H 2 / 2 ]T -^=(|n0) o6 + |0n) ob ). (9) 



It is straightforward to realize that the incommensurate 
nature of the Rabi frequencies Q(x(t), 0)y/n at which en- 
tanglement is transferred to each subspace spanned by 
{\n — 1, 1) , \n, 0)}ki,j generates quantum interference ef- 
fects. The result is that the amount of entanglement 
effectively passed from the ECS resource to the qubit re- 
ceivers is less than what is achieved using an entangled 



single-photon state. In order to see this, we have quanti- 
fied the entanglement between the flying qubits (initially 
prepared in their ground state), after their interaction 
with the driven cavities. The reduced two-bit density 
matrix is 



P12W 



/o 

A B C 

B A C 

\0 C* C* \-2A) 



(10) 



with 



„2n " 



CV — 

A = til E — E (C n m ) 2 R m T n - m sm 2 [ti(t)Vr 



n=l 



m=0 



B=NlTa 2 sin 2 [/iWL 



G 



iN 2 a {2Vfa^{t)] +EE 



x sin[/i(i)v s — m + 1] cos[/x(t) yjs — m]} 



and N a — A/" Q e~' a ! 2 / 2 , C 



£^2s+l j^Twjps—m-^-l /2 

~i ^0 m! V( s + 1 - m ) ! ( s - to ) ! 

(11) 



y n\/\m\(n — m)!]. 
In order to quantify the qubit entanglement, in the re- 
mainder of this paper we use the negativity of p\2 (t) [36[ , 



N[pi 2 {t)] =max[0,-2A_(*)], 



(12) 



where X-(t) is the smallest eigenvalue (with its sign) of 
the partially transposed two-qubit density matrix. 



A. Discussion of the results 

The behavior of the qubit entanglement under a 
NOON and an ECS driving and a perfectly unitary dy- 
namics is shown in Figs.[2ja) and[2fb) and[3]Ja) and[3jb), 
respectively. Clearly, a single-photon entangled state is 
able to prepare a sizably entangled qubit state depend- 
ing on the value of the frequency J7o- If the latter is 
such that the atoms perform an odd half-integer (inte- 
ger) number of Rabi floppings by the time they cross the 
field modes, the flying qubits leave the cavities with a sig- 
nificant degree of steady-state entanglement (in a separa- 
ble state), regardless of the motion (whether at constant 
velocity or accelerated) that they underwent within the 
resonators [cf. Figs. [2Ja) and^b)]. The working condi- 
tions of the process can be adjusted so that the values 
of the steady-state entanglement achieved in the two in- 
stances of atomic motion are comparable (cf. Fig. [2]). 
Needless to say, for free-falling qubits, a Rabi frequency 
much larger than the one that should be chosen under 
conditions of constant velocity is required, as the inter- 
action time between atoms and cavity fields is shortened. 
As anticipated, the use of an ECS driving does not result 
in an equally effective transfer efficiency. Figures [3ja) 
and Elb) show the negativity for both free-falling and 
constant-velocity qubits, optimized numerically over the 
parameters of the interaction and the amplitude of the 
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FIG. 2: (Color online) (a) Entanglement transfer for free- falling atoms and a single-photon entangled driving field (i.e., NOON 
states with N = 1), x(t) — xo+gt 2 /2, where xq = — 4o;o, <? = 9.82 m/s 2 , and for fio = 4 kHz (dashed line); S7q = 5.9 kHz 
(solid line), (b) Entanglement of atomic qubits crossing the cavities under conditions of uniform motion [x(t) = xo+Vt with 
xo = — 4ojo and V = 0.001 m/s] and with the same driving field used in (a). We have taken fio = 155 Hz (dashed line) and 
flo = 365 Hz (solid line), (c) Comparison between the unitary (solid line) and dissipation-affected (dashed line) entanglement- 
transfer performance for free-falling motion of the atomic qubits with Qq = 5.9 kHz. In the dissipation-affected simulation, we 
have taken fio/F ~ 10 2 . In all our calculations we have taken the cavity mode ^2,0(2;, 0) (i.e., the TEM20 mode), luo = 10 fim, 
and T = 0.9, and the qubits are prepared in |00) 12 ■ 
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FIG. 3: (Color online) (a) Entanglement transfer for free-falling atoms and ECS external driving with xo = — 4cjo, a=l.l, and 
flo = 6.1 kHz. (b) Entanglement for atomic qubits crossing the cavities under uniform motion (xo = — 4loq, V = 0.005 m/s) and 
flo — 850 Hz, under the same driving as in (a). The atomic steady-state entanglement at the exit of the cavities for both motions 
is comparable, (c) Comparison between the unitary (solid line) and dissipation-affected (dashed line) entanglement-transfer 
behavior for uniform motion of the atomic qubits (V = 0.005 m/s), Qo = 340 Hz, and a = 0.9. We have taken Qo/T ~ 10 2 . 
For all plots, we take the cavity mode ^2,0(2;, 0), too = 30 /jm, T = 0.9, and initial atomic ground state. 



ECS state (at a~l.l the entanglement between the flying 
qubits is maximum, uniformly with respect to the choice 
of the other parameters of the model). 



As seen from the analysis above, interaction times of 
the order of 10" 3 (for accelerated qubits) or 1CT 2 (for 
atoms moving at constant velocity) are needed in order 
to cross the cavities. These values of the interaction time 
might put the dynamics of the system well within the 
typical time scale of leakage of the cavity fields due to 
finite cavity quality factors. Therefore, in order to make 
our analysis consistent with a realistic experimental sit- 
uation, we should abandon the perfectly unitary descrip- 
tion adopted so far in favor of an open-system dynamics 
that includes the effect of field dissipation. This is done, 
in what follows, by implementing a quantum Monte Carlo 
unraveling of the system's dynamics [37j implemented by 
letting the system evolve through the non-Hermitian op- 



erator 

2 

H(t) = J2 H jkj (*) - iTktkj {h =A,k 2 = B) (13) 
j'=i 

with r the damping rate of each field (assumed to be 
the same for both the cavities, an assumption that can 
be easily relaxed if needed) and considering the effects of 
the quantum jump operators vTfcj-. The results of our 
numerical simulations of a large number of possible dy- 
namical histories of the system are presented in Figs.EJc) 
and[3Jc) for max t Q,(x(t),0)/T ~ 10 2 , which is a regime 
that can be achieved experimentally [2(| H|[ [in Fig. [3jc) 
we address the suboptimal case of a — 0.9 simply for 
convenience of calculations in light of the size of the trun- 
cated Hilbert space within which we have performed our 
quantum-jump calculations]. The process shows a notice- 
able robustness to the effects of field damping, leaving the 
qubits at the output of the cavities well entangled, irre- 
spective of the details of the atomic motion and the form 
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of the driving field. 



IV. IMPOSSIBILITY OF ENTANGLEMENT 
TRANSFER IN DISSIPATION-DOMINATED 
REGIMES 

Here we address the case of a dissipation-dominated 
dynamics of the cavity fields to demonstrate that no en- 
tanglement can be transferred in this case, regardless of 
the flying or trapped nature of the atomic qubits. More- 
over, besides determining the regime for effective entan- 
glement transfer, this approach will help us understand- 
ing some of the results that have been gathered in the 
previous sections. Let us start by introducing damping 
of the cavity fields into the dynamical evolution of the 
system, which will now be ruled by the master equation 

QtPVXAB = ~i[ J] Hjkj ) Pl2Ag] + ^ £j{p\2AB) 
i=l,2 j=A,B 

= (C + C)(pi2ab) 

with Pyiab the density matrix of the whole cavity-qubit 
system and (j = A,B) 



C j (a)=T(2jaj1 



(15) 



which is the Liouvillian describing the damping of each 
cavity field due to a low-temperature bath (as is typ- 
ical at optical frequencies) and acts on a generic den- 
sity matrix a. We now assume T 3> |f2(x, i.e., the 
bad-cavity limit, and trace out the cavity fields, so as 
to find an effective dissipative dynamical map for the 
qubits only, driven by a structured quantum environment 
that exhibits quantum correlations. Such an effective 
description is gathered following standard strategies for 
the derivation of adiabatically eliminated master equa- 
tions 139] as 



d t pi2 = Tr A s < C 



e Ct C (p s I B 



pi 2 )dt 



(16) 



which has been obtained by defining a projection oper- 
ator V such that V p\iAB — Pab ® /°i2 with p" B the 
steady state of the cavity fields and p\ 2 the reduced den- 
sity matrix of the atoms only. The projection opera- 
tor is such that VC^PpxiAB = 0, which is a key prop- 
erty for the derivation of Eq. (p~6|) . In Ref. [12[, it has 
been shown that the explicit form of the reduced mas- 
ter equation depends on a Kossakowski matrix that is 
fully determined just by the second moments of the en- 
vironmental modes. We introduce the covariance matrix 
M(pab) with elements Mi j (pab) = Tr[pAB{4i, Qj}/2], 
where q — (xaPaxb Pb) is the row vector of the two- 
mode field quadratures. Any two-mode covariance ma- 
trix can be transformed, by means of local symplectic 
operations, into the form [401 ] 



, n c 
M = \ T 

, c rn 



(17) 



where p = pl 2 (p = n, m) account for the variances of 
each local mode and c = diag[ci,C2] describes the inter- 
mode correlations. In our formalism, the two-mode vac- 
uum state corresponds to n = m = 1/2 with c\, 2 = 0. 
For a generic covariance matrix M, the Kossakowski ma- 
trix that fully determines the dynamics of the two-atom 
system can be constructed as 



K = 7 



n + m 



m + iH 



(18) 



with 7 cx \Cl(x, y) | 2 /r the effective two-qubit coupling 



rate and £ = 



1 

-1 



the single-mode symplectic ma- 



trix [12J, |16[ . The reduced dynamics of the atoms is thus 
built as 



dtP\2 



4 



K aP (O a p 12 O p - {O p O a , Pl2 }/2), (19) 



'a-2 



for 



where O a = a\ (8> 1 for a = 1 , 2 and O a = 1 < 
a = 3,4. Here, o"^^ is the x(y)-Pauli matrix of qubit 
j = 1,2. Complete positivity of the map originating 
from Eq. (JTHJ) is ensured for K > 0, which is equiva- 
lent to the Heisenberg-Robertson uncertainty principle 
for the covariance matrix M [3{|. In what follows, we 
will consider both the temporally-resolved two-atom dy- 
namics achieved by solving Eq. (|19p and the steady-state 
one obtained by setting dtpi 2 = 0. 

We start with the NOON driving field. Upon ex- 
plicit calculation, it is straightforward to check that 
MQonoon) = (iV-h 1)1/2 for any TV > 1 (here, pnoon = 
|NOON) (NOON|). This simply implies that all the cor- 
relations in a NOON state with more than one photon are 
encoded in higher-order moments of the mode quadra- 
tures. As the process of entanglement transfer ruled by 
the interaction Hamiltonian H\a + H 2 b relies on the sec- 
ond moments only, a point that has been duly stressed 
in Ref. [13], it is clear that such a driving field does not 
embody a useful resource for the process, regardless of 
the dependence of the coupling strength £l(x, y) on the 
details of the atomic transition across the cavities. In 
turn, this provides a rigorous explanation of the obser- 
vations made earlier for the case of unitary dynamics. 
For N = 1, on the other hand, the covariance matrix is 
nondiagonal and reads 



M(piooi) = - 



/2 1 0\ 

2 1 

10 2 

Vo 1 2/ 



(20) 



This gives rise to the steady-state two-atom density ma- 
trix 



P ss = tt: 



1 

12 



'1 ON 

2 10 

12 

v0 7/ 



(21) 
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FIG. 4: (Color online) (a) Temporal behavior of the elements of the two-qubit density matrix for the case of an entangled 
single- photon driving field (solid lines). The dashed lines are the corresponding steady-state values. We have used the notation 
(pi2)ij,ki = (ij\pi2ss\kl) . The initial state of the atoms is taken to be |00) 12 . (b) Quantum discord versus dimensionless 
interaction time yt for the case addressed in (a), (c) Atomic steady-state discord against the amplitude a £ K in an ECS 
driving the entanglement-transfer process. Dynamically, the two-atom state is always separable with discord following a trend 
(at a set value of a) very similar to what is shown in (b) for an entangled single-photon driving field. 



which has positive partial transpose, as it is straightfor- 
ward to check, and thus describes a separable state. In 
Fig. BJa) we plot the time behavior towards steady state 
of the elements of the two-atom density matrix. Interest- 
ingly enough, not only is the two-atom state fully separa- 
ble at all instants of time of the dynamics, it is also very 
weakly quantum correlated, in general. In order to sup- 
port our claim, we have computed the entropic version 
of quantum discord [4l[ proposed in Ref. [42j to test for 
more general quantum correlations seeded into the state 
of the two atoms. Quantum discord strives at capturing 
quantum correlations of a broad nature in a multipar- 



tite quantum system, witnessing the non-classicality of 
the way correlations are shared by the elements of (in 
general) a many-party register. In our case, a non-zero 
value of quantum discord in the bipartite atomic system 
at hand signals the existence of quantum features in the 
density matrix describing the state of atoms 1 and 2. The 
results are shown in Fig. [4fb) , where it is shown that a 
very small degree of discord is shared by the atoms. The 
steady state, though, remains non-classically correlated, 
although separable. 

We now consider an ECS driving field [cf. Eq. (j6|)], 
whose covariance matrix is 



m(pecs) = M a 



( a 2 + (l + e- a2 )(a 2 + 1) 

e- a \l-a 2 ) + l 



e~ a a 2 a 2 + (1 + e- Q ")(a 2 + 1) 



V 







e~ a a 2 



(22) 



a (l-a 2 ) + l J 



While both steady-state and time-resolved density ma- 
trices can be computed analytically using the formal ap- 
paratus described above, their expressions in terms of 
the entries of the covariance matrix in Eq. (|2"2"j) are too 
lengthy to be reported here. We can then check for insep- 
arability of the two-atom state, finding that, as in the en- 
tangled single-photon case, no entanglement is set by the 
dynamics, regardless of the form of the Rabi frequency 
£l(x, y), both dynamically and at the steady state. As for 
the quantum discord, we have computed it for the steady- 
state density matrix against the amplitude a entering the 
driving field [see Fig. H[c)]: the discord is nonmonotonic 
against the ECS amplitude, showing that (very modest) 
non-classical correlations are set preferentially at small 
values of a (the maximum of the discord being found for 
a ~ 0.64). As a increases, the discord vanishes, thus 
giving back a state with only classical correlations, at 



most. In fact, for a 3> 1 the two-atom state is com- 
pletely uncorrelated. Indeed, the off-diagonal terms in 
the steady-state density matrix go to zero as a increases, 
thus leaving a diagonal state with no correlations at all 
(actually, the atomic state tends towards a maximally 
mixed one). 

The reason for the trends highlighted in this section is 
very clearly related to the fact that the covariance ma- 
trix of the field states driving the cavities fails to vio- 
late the criterion of positivity of partial transposition. 
Needless to say, as the field resources are non-Gaussian, 
this implies only that correlations are encoded in higher- 
order moments of the field's quadratures, as mentioned 
above. The dissipation-led entanglement-transfer pro- 
cess, though, strongly relies only on the mentioned second 
moments: as far as the process at hand is concerned, us- 
ing the driving fields analyzed above is not different from 
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driving the cavities with separable Gaussian states. This 
analysis suggests the pathway that should be pursued, 
experimentally, when entangled single-photon and ECSs 
are used as resources: an almost unitary process needs to 
be in place, in order for the transfer process to actually 
take place. In contrast, should the transit of the atom 
occur in times comparable to the cavity lifetime (alter- 
natively, should the quality of the cavities be not high 
enough), the scheme will fail, leaving, at its best, atomic 
states that are only weakly quantum correlated. 

V. CONCLUSIONS 

We have investigated the entanglement transfer be- 
tween flying qubits and optical fields based on a cavity- 
QED system. We have considered two typical quantum- 
correlated driving fields, NOON and entangled coherent 
states, and found that entangled atomic qubits can be 
prepared, effectively, regardless of the details of their mo- 
tion. We have discussed diverse aspects of our system, 
including its robustness to the effects of non-negligible 
cavity dissipation and sensitivity to the value of the Rabi 
frequency, finding that the maximum entanglement the 
qubits can achieve is not determined by their motion 
but by the characteristics of the cavity fields, which in 



turn depend on the external driving and cavity damp- 
ing rate. The high degree of robustness of the scheme 
(for experimentally realistic parameters) allows consider- 
ation of the consecutive passage of pairs of independent 
qubits through the entangled cavities and the construc- 
tion of a stream of entangled qubits. The potential and 
limitations of such scheme will be addressed elsewhere. 
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